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We present a scaling technique which transforms the evolution problem for a non- 
linear wave equation with small initial data to a linear wave equation with a distribu- 
tional source. The exact solution of the latter uniformly approximates the late-time 
behavior of solutions of the nonlinear problem in timelike and null directions. 
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<C ' I. INTRODUCTION 

t> 

For quite some time it is known that small global solutions of nonlinear wave equations 

ic£ Du = F(u) (1) 

develop late-time power-law tails. Several authors [l|, 0, H, 0] have given pointwise decay 
estimates which imply bounds on the decay rates as t — > oo. However, they do not show 
that these decay rates are optimal. It was only recently that we were able to calculate the 
exact decay rates and gave conditions on the initial data under which the decay estimates 
are optimal 0, lEf. Our proof was, however, restricted to the asymptotics in the timelike 
^ . direction, i.e. t — > oo with x being fixed. Here, we give a uniform asymptotic formula which 
OO '. holds in the timelike and null infinity, i.e. for t — > oo and < \x\/t < 1. 

We consider the initial value problem for nonlinear wave equations of the form (pQ) with 
small initial data 

O 

in three spatial dimensions, i.e. u : (t, x) G M+ x R 3 — ► R, where e is a small number. The 
nonlinearity behaves like 



{u,d t u)\ t =o = (ef,eg) (2) 



X ' F(u) = u p + 0(u p+k ) for u -> (3) 

with p > 1 + y2 and k > (if p is non-integer the expression u v is to be understood 
as |m| p_1 m everywhere where needed). The main idea developed in this article is based on 
the observation that in the limit e — > the solutions u £ tend under suitable scaling to 
some nontrivial u* which satisfies a linear wave equation with a distributional source. This 
equation can be solved exactly. We show that for small but finite values of e the solutions 
u e are near to u*, in a suitable sense with a uniform error bound, such that u* determines 
the late-time asymptotics of u £ and plays a role of a universal attractor in the evolution. 
Our main result is the asymptotic formula 



u(t, r) = euo(t, r) + e p — 

r 



(t - r)P~ 2 (t 



\p-2 



+ O p _x {e p+x ) (4) 



where Uq is a solution of the linear problem with data d2J) and the second term represents 
the leading nonlinear asymptotics with A p being the only trace of the initial data. The error 
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term O p _i (to be made precise below) is uniformly small relative to the leading asymptotics 
as e — > and A is some positive number depending on p and k. 

This asymptotic formula recovers the known timelike asymptotics u ~ e p '/t p ~ 1 |5j as well 
as null asymptotics u ~ 1/t and is regular at r = 0. 

The idea of scaling has been inherited from the Doctoral Thesis of Hans Lindblad (fif where 
he, in contrast, used similar scaling technique to treat blowup in finite time of solutions to 
the wave equations Hu = u p with p < 1 + \/2. 

The scaling method for small global solutions introduced here can also be straightfor- 
wardly extended to more complicated systems of nonlinear wave equations as long as appro- 
priate a priori decay estimates are known. However, even if no such estimates are available 
this technique can act as a (nonrigorous) guide to the late-time asymptotics of the studied 
system. This method can be also used to study the sub-leading asymptotics, order by order 
in an iterative way, which we plan to describe in a subsequent publication. 

For more transparency in this introductory paper we restrict ourselves to spherical sym- 
metry as being known to dominate the late time asymptotics of the studied equation anyway. 
For generalization to full 3-dimensions, which here is a purely technical element, one can 
repeat the corresponding steps from Lindblad 0]. 

This paper is organized as follows. First we introduce the idea of scaling and discuss the 
limit e — > in a non-rigorous way. Then we prove a theorem on late-time asymptotics of ([1]) 
by comparing the original system with the limiting rescaled one. Finally, we discuss some 
applications as well as the extensions of the technique. 



Notation 

With the symbol (x) :— 1 + \x\ we define a space-time weighted-L°° norm 

IMU~ : = IK* + \x\)(t~ \x\) p - l u{t,x) \\ L oo (M Wy (5) 

Its finiteness guarantees the decay of u like 1/t on the lightcone t ~ \x\ and like l/t p for 
fixed x as well as for fixed t. We also define a smaller norm restricted to the region 

Q := {(t,x) G M^+ 3 :t- \x\ > 1} 

:= IK* + NX* - M) p ~Vt,x)|| L oc (n) < |H| L oo. (6) 

In 3 spatial dimensions by IR^ +3 we will understand R + x M 3 and in spherical symmetry 

will mean M + x R + with {0} G R + . 
By D -1 we will denote the inverse of the wave operator □ having the properties 

u = n~ 1 F Uu = F and (u, d t u) | t=0 = (0, 0) (7) 

In 3 spatial dimensions is a positive measure what gives a useful comparison property 

f > g n-*F > D^G. (8) 

We define the symbol O q as an extension of the big-0 symbol for e —>■ to the space L£°, 
i.e. for a family F e e and G G C(R), g{0) = we say 

F £ = O q {G{e)) :^ \\F £ \\ LT =0{G{e)). (9) 

Whenever we write \3~ 1 0(F £ ) or n~ 1 O p (s q ) we mean D _1 is acting on a function of pre- 
scribed decay which is at least in C 2 (R. + +3 ). 
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II. THE IDEA OF SCALING 

In this section we want to sketch our strategy of using scaling to obtain exact late time 
asymptotics for small data. To this end we consider a simplified version of (CQ) 

Uu = u p (10) 

with integer p > 3 and initial data 

u(0,r) = ef(r), d t u(0, r) = eg(r) (11) 

being smooth functions of compact support in [0, R\. First, we solve the linear equation 
with removed scale factor e 

Du = 0, (12) 

«o(0,r) = /(r), dtUo(p,r) = g(r). (13) 

Its solution can be written in the form 

. . hit - r) - hit + r) 

u {t,r) = ^ >- ^ J - (14) 

r 

where 

x 1 f°° 

K x ) : = 'of^ ~o y9(y) d v ( 15 ) 

has support in [—R,R] (the functions f(r),g(r) have been symmetrically continued to neg- 
ative r). Next, we subtract the linear solution from the nonlinear by introducing 

w(t,r) :— u(t,r) — suo(t,r) (16) 

which satisfies 

Uw = (w + eu ) p . (17) 

Now, we scale this function to 

W E (t, r) := e~ b w{e- a t, e~ a r) (18) 

with b = p + a(p — 1) and some a > to be chosen later. It satisfies 

DW £ (t,r) = e- a [e {p - 1)+a{p - 2) W £ (t,r) + e- a u (e~ a t,e- a r)] P 
= e~ a [e- a u {e~ a t,e- a r)] p 

- a r)] p - k e^-V +a( *-^ k W e k (t 7 r). 



+ £ ~ a [^ a Me~Xe- 

k=l ^ ' 



(19) 



For this equation we want to consider the limit e — > 0. From our previous work [7] we have 
global weighted bounds on u and w when e is sufficiently small 

Ce Ce p 
^)\< {t + r){t _ r)p _ 2 , Ht,r)\< {t + r){t _ r ^ r (20) 
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It implies that the functions u, w and by (fl8"|) also W £ (t, r) are bounded in the weighted-L 00 
norm 

\H\l~, < Ce, I^IU^ < Ce\ < Ce'^, \\W e \\ z?Li < C. (21) 

Let us recall the following fact of the distributional calculus: any smooth function H(x) 
of compact support can be squeezed to the delta distribution under appropriate scaling as 
e -> 0, 

e- x H{e- x x) -> C H 5(a;), C H := I H(x) dx. (22) 



Having this in mind we observe that all terms in ([19]) containing powers of uo, by use of the 
representation ffl4|) . will have a distributional limit 

e~> [e-'Me-Xe-'r)]'-* c j^ ~ + r) (23) 

with 

C fc : = / Z^cfe. (24) 



The term 5(t + r) will further play no role since its support is outside of the region of our 
interest t + r > 0. Applying these limits to equation (fl9l) we observe that the first term 
tends to 

while the sum over k vanishes as e — > because it converges to the delta distributions times 
functions W £ (t,r) which are uniformly bounded by Ce~ a ^ p ~ 1 ^ k and factors e[(p~ 1 )+ a (p~ 2 )] fc ) 
combined together giving 

£ [(p-l)(l-a)+a(p-2)]k _^ q (26) 

for any p > 2 and < a < p — 1. This leads us to a limiting equation 

DW(t,r) = C p 6 -^l (27) 



which can be solved exactly 

9(t-r) 



W(t, r) = A 



p 



r L(t-r)P~ 2 (i + r)?- 



(28) 



with A p := 2 p ~ 3 C p /(p - 2). After W is traced back to the original solution u it will provide 
the leading asymptotics for small e. In the rigorous treatment all error terms which we 
ignored here must be estimated as uniformly small is some region of spacetime. The crucial 
point will be to show that W e tends to W in the norm L^_ x , i.e. 

\im\\W £ -W\\ Too =0, (29) 

hence 

W £ (t,r)=W{t,r) + O p _ 1 (l) (30) 
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uniformly in t — r > 1 . Then we can reconstruct the asymptotics of u 



u (t, r) = eu Q (t, r) + e p+a{p - 1] W £ {e% e a r) 



eu (t,r) + e p A, 



Q(t-r) 



(t-r) 



P-2 



(t 



\p-2 



O p _i(e* 



(31) 



which holds uniformly (in the sense of the L^^-norm) in t — r > e~ a . 

This asymptotic formula is richer than the one we obtained in [5| for it additionally gives 
an analytic dependence on r in the whole region t — r > s~ a . Hence, it can be thought of as 
asymptotics for t — > oo including dependence on the parameter r/t in the range [0, 1). The 
asymptotic formula of [5[ is reproduced by fixing r and letting t — > oo what gives 



u(t,r) ~ e p /t 



p-i 



(32) 

r and vanishes for t^> r 



because u is localized in the vicinity of the outgoing lightcone t 
leaving only the contribution of W. 

Here, essential is that the error term represented by O p -\(e p ) is controlled in the L^_ 1 
norm, because only this guarantees that the error has at least the same falloff as W while 
having a smaller amplitude for e — > 0. 



III. MAIN THEOREM 

In this section we proceed with allowing general nonlinearities of the form ([3]) with any 
real p > 1 + y/2 (not necessarily integer). The initial data will still be of compact support 
[0, R], although data of sufficiently fast fall-off at infinity could be also considered within 
this scheme. This would, however, lead to somewhat more complicated asymptotics as an 
incoming radiation may alter the solution at late times (for more details on competition of 
the asymptotic effects, see 0]). We choose (/, g) G C 3 x C 2 such that the following norms 
are finite 

/o := IK^rVWIUoc, h := \\(r) p Vf(r)\\ L ~, g := \\(r) p g(r)\ L ~. (33) 



By the results of [4j], it guarantees existence of a classical solution u G C 2 (M 1 f +1 ) fl L™ x for 
small values of e which multiplies the initial data in (J2J). The regularity can be relaxed to 
weak solutions in the energy space u G C°(1R, H l ) fl C 1 (1R, L 2 ), but the estimates will then 
become much more delicate and the pointwise asymptotics only true almost everywhere (cf. 
[3] for such analysis). 

Since we are interested in the late-time behavior of the solutions we must prescribe the 
asymptotic region in which our approximation should hold as e — > 0. For there is some 
freedom in doing this a free scaling parameter a will be introduced. The quality of the error 
bound will then depend on its value, i.e. on the rate the asymptotic region shrinks to the 
vicinity of infinity when e — *■ 0. 

Theorem 1. The solution u of the initial problem (ITT) - ([21 has the following asymptotic 
behavior for e — > and a given scaling parameter < a < p(p — l)/(p + 1) 



^4 

u(t, r) = eu (t, r) + e p — 

r 



(t - r)P- 2 (t 



\p-2 



+ O p _x (e p+x ) (34) 
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where Uq solves the linear problem (|T2"1) -(IT3" | ). A p := 2 P ~ 3 C P / (p — 2) with C p defined via ((24 
and ( TT5l) and A := min{[p(p — 1)(1 — a) + a((p — l) 2 — 2)]/p, k(l + a), a}. TTie error term 
is defined in (J9j) and means that the asymptotics holds w.r.t. to the weighted-L p x l 1 norm. 
Here, it is restricted to the region t — r > l/e a what implies a uniform convergence there. 



Before we prove the Theorem, we first cite a useful lemma 
Lemma 1. If p > 1 then for any real numbers a, b 

\a + b\ p < 2 p - 1 (|a| p + \b\ p ) 

and 



a + b) p -a p -W < 2 p ~ 1 p(\a\ \b\ p ^ + \b\ \a 



p-i 



(35) 



(36) 



The inequality remains true when all expressions of the type x p on the left-hand side are 
replaced by \x\ p or |a;| p_1 x. 

Proof. See e.g. Lindblad 0, Lemma 8.3] with minor modifications not affecting the result. 

□ 

Proof of TheoremQl From [4, Th.2] we know that there exists a classical solution u 6 
C 2 (IR+ +1 ) which belongs to the L^L 1 space. 

Let Uq solve the linear equation ffT2|) with the initial data ([TBI satisfying (133]) . Then, by 
0, Lem. 4a], there exists a classical solution u G C 2 (1R 1 I +1 ) which satisfies 

IMIx-, <C-(g + fi + fo). (37) 

(This can be also directly verified by using the explicit form (fT4"j) of Uq.) Define 

w(t,r) := u(t,r) — eu (t,r) (38) 

and 

W e {t, r) := e- b w(e~ a t, £" a r) (39) 
with b = p + a(p — 1). From Theorem 3.1 and inequality (3.7) of [7] it follows 1 that 

\\u\\ L?Li <Ce, \\w\\ L?Li <CeP. (40) 

It implies 

<Ce~ a ^ (41) 

and 



\W f 



ellL. 



p-i 



£ -[ P+ a( P -l)] || {t + r) {t _ r y-2 w{ t/ £ a^ r/e «) \\ Loo(t _ r>1) 

(t + r)(t- r) p ~ 2 ■ (t/e a + r/e a )(t/e a - r/e a ) p ' 2 \w(t/e a J r/e a )\ 



e 



-b+a(p-l)] 



-[p+a(p-l)]| 



p-l 



(t/e a + r/e a ) (t/e a - r/e a ) p ~ 2 

£ «(P-1)(1 + t + r )(l + t _ r )P-2 



L°°(t-; 



(e a + t + r)(e a + t - r) p ~ 2 



L°°(t-r>l) 



p-i 



2P-l £ -a(p-l) 



< 2 p - 1 C, 



(42) 



1 Notice a slightly different notation: uq here corresponds to u\ in 
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where we have used the inequality l + t±r<2(e a + t±r) which holds for t — r > 1. 
W £ satisfies the following nonlinear wave equation 

n\V e {t,r) = e -P-«0H-i)_p ( £ v+a{p - l) W £ (t,r) + Eu (e- a t,e- a r)) 



(43) 



Now, lead by the observations made in the previous section, we introduce the limiting 
function W(t,r) defined by the limiting equation 



C k ■ = / h k (x)dx 



where 



and show that indeed W £ converges to W in norm Li x L 1 , i.e. 



lim||W s -H^||foo =0. 



(44) 
(45) 

(46) 



Note, that it cannot be shown that the right-hand side of fj43|) converges in L^-norm to the 
right-hand side of (|44p as the norm of the latter is infinite. Convergence of W £ to W, after 
inverting □, does not work either, since the L^-norm of W is infinite. However, both W £ 

and W have finite L^ c L 1 -norms. This reflects the fact that pointwise convergence happens in 
the region t — r > 1 but not in the vicinity of the light-cone t ~ r. 

Let us consider their difference and split it into four parts for further analysis 

\W £ -W\(t,r) = p-i £ -P-»(p+i) F ( £ P+ a ^W £ (t,r) + eu (e- a t,E~ a r)) - n~V| 
< e- a D- 1 (e {p - 1)+a{p - 2) W £ {t,r) + e- a u {e- a t,e- a r)y 

e [{p - 1)+a ( p - 2)] W £ (t,r)Y - (e- a u (e~ a t,e- a r 



+ e 



-p-a(p+l) 



p+k 



(47) 



□- x o( [e p+a{p - l) W £ {t,r) +eu (e- a t } e- a i 

£ p(p-l)+a[(p-l) 2 -2] | jyP(£ ; r ) | 

+ |n- 1 e-° (p+1) «g(e- a t,e-°r) - D^VI 
=: A p + A + A w + A M . 

The first term, A p , will be first estimated algebraically with Lemma [1] and then in another 
Lemma, by Holder inequality, reduced to a product of two separate expressions containing 
□^H 7 ] 3 and D -1 ^. Then two next Lemmas will transform these terms to norms of W F 



and uq. Finally, use will be made of the norm bounds (|4T|) and (!3Tj) . The second term, 
Ao will be estimated similarly. The third term, A\y, will be estimated by the same means, 
using norm bound (j4T|) only. Finally, the estimate of the last term, A M , which is the core of 
the technique, will be based on the fact that Uq scales to a distribution, as observed in the 
previous section (cf. fT23"j) ). 
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Lemma [T] together with the positivity of □ 1 (cf. (JSJ)) give 

A p < e^U- 1 (e {p - 1)+a ( p - 2) W £ (t,r) + e~ a u (e- a t,e- a r)y 

_ r £ i(p-i)+a(p-2)] w ^ r j\ p _ f £ - a Uo (e- a t,e- a r) 

< P 2 p ~ 1 e- a D- 1 ( e Kp - 1)+a{p - 2)] W e (t, r) e- a uo(e- a t, e' a r) 



p-i 



+ 



p-1 



P 2 p - 1 [e p - 1 - 2a D- 1 ( W £ (t,r 

_|_ e (p-l) 2 +ap(p-3)Q-l 



p-i 



e- a u (e- a t,e~ a r) 
p-i- 



) 



(48) 



W e (t,r) e- a u (t,s- a r) 



Now, we are going to separate the products under □ . 
Lemma 2. For positive a, /3,p > such that a + (3 = 1 we have 



(49) 



Proof. The proof is a straightforward calculation based on the fact that D -1 / can be written 
as a convolution E * f, where E is the fundamental solution of □, and use of the Holder 
inequality with powers 1/a, 1/(3 (cf. again Lindblad jl, Lemma 8.11]). □ 

We apply this Lemma to our estimate on A p above with / = |W e |,(7 = |w | an d ct = 
l/p,f3 = (p — l)/p or vice versa. Then we get 



A p < p2 p 



-i 



£ p-i-^( a -i\ W p (t, 



i/v 



□- 1 |ug(e-%£-°r) 



(p-i)/p 



_|_ £ (p-l) 2 +ap(p-3) ( 



(p-l)/p 



!<(£"% e~ a r) 



i/p-i 



(50) 



It is convenient now to switch from pointwise to norm estimates. We get 



\W £ -W\\ Tac < ||AJ| £o o 

p— 1 p— 1 



|A || foo + ||A 



+ IIA 



(51) 



p— 1 p— 1 

These weighted-!/? 3 norms have the following property (cf. Asakura [1, Eq. 2.32] 

Lemma 3. Let v,w G for some q > 1. T/ien /or any positive a,j3 > snc/i £/ia£ 
a + /3 = 1 we have 



\v a w l3 \\ L oo < \\v\ 



\w 



1/3 



T/ie same /io/ds for the L°° norms. 

By this Lemma with q = p — 1 we get 



p-i 



< p2 p - x 



p-l-2o 



1/p 

2s 



(p-l) 2 +ap(p-3) 



'P-1 

□ _1 lw? 



□- 1 |uS(£- t,e- r) 



(p-i)/p 



(p-i)/p 



□-^(e-^e^r)! 



i/p 



>-i 



(52) 



At this stage we need two estimates for the two different terms containing W e and uq, 
respectively. The first we cite from Asakura [1], Cor. 2.4 and Eq. 2.33] 
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Lemma 4. Let v G C 



2/Dl+ll 



fl L£° for some q > 1 . TTien /or any p > 1 + y/2 



|p-V|P|| LOO < C\\v\\ p 

II I I ll^g II l\L, Q 

with some C > provided q < p — 1 . 

It holds also when is replaced by on the left-hand side (but not on the right-hand 
side!). The proof is almost identical. 

The second Lemma, using explicitly the scaling properties, we prove below 

Lemma 5. For Uq given by (|14p and for any q > 1 



□ _1 \u q {e- a t,£- a r)\ 



(53) 



J q-1 



with some C > 0. 

Proof. Observe that due to scaling properties of the operator D" 1 we have 
[□"^(e-., £""-)] (t, r) = e 2a [D" 1 ^] {e~% e -r). 



(54) 



From Lemma H] and the bound (j3~Tj) we get ||D 1 |wo| 9 ||£°° < C|l M o||z / °o i < C- It gives the 
pointwise estimate 

/ x C 

(t,r) < 



for t — r > 1 which scales to 



e- a t,£- a r) < 



(t + r){t-r)i- 2 
Ce a( q -i) 



< 



C ' £ a(q-1) 



(e a + t + r)(e a + \t - r\)i- 2 ~ (t + r) (t - r)i~ 2 
for t ± r > 1. Combining (1541) with fl56l) we get the thesis. 
Now, Lemma H] and the bound (jUJ) give 



< \\W 6 \\ P TOO < Ce- a{p - 1)p 



while Lemma [5] gives 



kig(£- a t,£- a r) 



Finally, we arrive at 



< C 



p-i 



(p-l)(l-a)+a[(p-l) 2 -2]/p + £ (p-l){(p-l)(l-a)+a[(p-l) 2 -2]/p} 



(55) 

(56) 
□ 

(57) 
(58) 

(59) 



with some constant C > 0. Since we assumed p > 1 + y/2 it holds (p — l) 2 — 2 > and both 
e's are raised to a positive power (> y/2) for any < a < p(p — + 1). The power of 
the second term is always bigger than that of the first term for their ratio is p — 1 > 1. 
For the Ao term we again use positivity of D _1 and Lemma [T] to get 



|A | < 2 p+fc - 1 n- 1 O p _i^ p( ^ 1)+a[( ^ 1)2 " 2l+[p+a(p " 1)]fc |W / e (t,r)| p+/c 



+ e- a(p+1)+k \u (e- a t,e- a r 



a r \\P+k 



(60) 



10 



and after taking the norm 

||A ||roo = o ( £ p(p-VM(p-v 2 -2}+\p+»(p-m\\n- 1 \w £ {t 

+ O (V«(p+!)+ fe || Q-i \u {e- a t, e- a r) \ p+k \\ Zoo i 
Using the results obtained above we get 



\p+k\ 



(61) 



IA 



o || £oo _ Q ^ £ p(p-l)(l-a)+a[(p-l) 2 -2]+kp^ + q ^ £ fc(l+a)J _ 



(62) 



It is now straightforward to estimate the Ayy term. By Lemma H] we have 



p-1 



< £ p( P -i)+a[( P -i) 2 -2] \\W e (t,r)\\ p L oc < Ce p{p ~ 1){l ~ a)+a[{p ~ l)2 ~ 2] 



(63) 



p-i 



where the power of e is again positive for any a > 0. 

Finally, we have to estimate the most important "error" term 

A M = \n- 1 s- a{p+1) u p {e- a t,e^ a r) - D^VI 



(64) 



which will tend to zero as e — > because the compactly supported outgoing wave Uq shrinks, 
under scaling, to the distribution \i localized on the lightcone t = r. First, observe that only 
the outgoing part of the radiation contained in uq contributes to the norm ||A 



a* II is 



since 



p- 1 e- a{p+1) u p {e- a t,e- a r) - CTV| 



n -i = - a He- a (t-r))-h(E- a (t + r))y _ ^ 



y-P 



< 



, y (E -(t-r))_ . 

r p ^ 



v p 



(65) 



and the last term (ingoing wave) vanishes when t — r > 1. The reason is that in the 
(spherically symmetric) Duhamel formula for \3~ 1 F with F(t, r) = e~ a h p (e~ a (t + r))/r p 



□ _1 .F(t,r) 



du 



t-r 



dv (u — v) 



e- a h p (e- a u) 
(u - u)f 



(66) 



the integration runs only over u > t — r and hence e~ a u > R is outside of the support of h. 

So we need to deal only with the first term (outgoing wave). We define an auxiliary 
function 

%{t,r,p) := D- 1 \h p (p) 6 t ! ~ ri 



r p 



(67) 
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and observe that 



! _ a h*(e-(t-r)) ! 
o„. := U £ U /i 



□ 



-i 



'^( e -V)- ( *" P " r) 



dp-e- a / □ 



i-i 



6{t - r) 



dp 



e a J $ p (t — p,r,e a p)dp-e a J $ p (t,r,e a p) <^P 

j dp j dr d T § p (T,r,e~ a p) 
J Jt-p 

-e- a [ dp h p (e- a p) [ dr 
J Jt-p 



[t — r 



[t + r 



where in the last step we differentiated the explicit expression for $ f 



$p(t,r, p) 



1 



1 



h p (p) 

p-2 [{t-r)P- 2 (t + r)P- 2 



Now, we can estimate 



/+e a R ft 
dp \W(s- a p)\ / dr 
■e a R Jt-p 



The inner integral can be bound 



dr 



t- P 
P 



(r - r)P- 1 
1 



(r + r)?- 1 
1 



< p sup 

t-p<r<t 



— r^P- 1 



[t — r 



[r + r 



(t - r - p) p - 1 (t + r - p)p~ 1 
C 

< Pit—. — i^t: — : — i^^tt < P 



<P 



[t -r)P- 1 
1 



[T + r 



\p-i 



1 



(t — r- e a R)P- 1 (t + r- e a R)p~ 1 



(t + r - e a R) (t-r- e a R)P~ 1 ~ (t + r)(t- r)P^ 



(68) 



(69) 



(70) 



(71) 



where we have used the facts that p G [— e a R, +e a R] and t ± r > 1 > e a R together with a 
simple algebraic inequality of Bernoulli's type (cf. proof of Lemma 1 in j8j). This leads to 



K\< 



e- a C 



(t + r) (t - ryp- 1 



-e a R 



e a R 



for t — r > 1 where 



Taking the norm gives 



L° 



dpp\W(e- a p)\ 



P l \h k (p)\dp- 



e a C'C Pt i 



(t + r) (t - r)P- 1 



(72) 



(73) 



(74) 



and we see that has faster decay than other error terms. In order to respect this fact 
we instead of using further the norm-estimate (I5T1) come back to the pointwise notation and 



12 



find 



\W. -W\< IAJ + |A | + |Aw| + |A 



r (p-l)(l-a)+a[(p-l) 2 -2]/p 



= O p _! (e 

+ o P (q}e°) . 

We see that all terms tend to zero for e — > when a > and < a < p(p — l)/(p + 1). This 
result can be written as 

W £ (t, r) = W(t, r) + O p _! (e Ao ) + O p (e a ) (76) 
as e — > uniformly in the region t — r > 1, where A > depends on the values of a and 



An 



k(l + a 

The exact solution of (144p reads 

6(t-r) 



p(p-l)(l-a)+q[(p-l) 2 -2] for P(p-l-fc) 
p ' p+l+pfc 

otherwise. 



(77) 



W(t,r) = A 



(t - r)P- 2 (t + r)P~ 



V(t,r)G 



>i+i 



(78) 



where A p := 2 P 3 C p /(p — 2). Finally, scaling all functions back, we find 



u(t, r) = eu (t, r) + e p+a{p - l) W £ {e a t, e a r) 

Q(t - r) 

= eu {t, r) + e p A p ^ '- 

r 



l(t-r] 



p-2 



(t + r. 



p-2 



(79) 



where the error bound is now uniform in the region t — r > e~ a . Observe that in that 
region we have O p (e p ) = O p _i (e p+a ) so the total error can be written as O p _i (e p+A ) where 
A := min{Ao, a} what finishes the proof of Theorem [TJ □ 



IV. DISCUSSION 



In the case when F(u) = u p the result is 



u(t, r) = euo(t, r) + e p A p 



Q(t-r) 



L(t-r) 



p-2 



(t + r)P" 2 _ 



+ O p ^(e p+Xo )+O p (e p ) (80) 



with A := [p(p - 1)(1 - a) + a((p- 1) 2 -2)] /p. The first error term describes corrections with 
the same decay in time as the leading order asymptotics but entering with higher powers of 
e while the second error term stays for a correction entering at the same nonlinear order (the 
same power of e) but having faster decay in time. By a suitable choice of a, and thus the 
asymptotic region of spacetime, one can make the first or the second error term dominant 
and explicitly evaluate its leading asymptotics from initial data. Hence, the technique can 
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be extended to study the sub-leading asymptotics order by order in the powers of e. By 
defining another auxiliary function v :— u — euq — e p w and scaling it properly and taking 
the limit e — > one can obtain another linear wave equation for v with a distributional 
source (containing now 5'(t — r)). This procedure can be, in principle, iteratively repeated 
to arbitrarily high orders in e, however the number of terms to be analyzed in the limit 
£ — > quickly increases. In a forthcoming publication we want to apply this generalization 
to the well-studied cubic wave equation and derive the late-time attractor found by Bizon 
and Zenginoglu in [9|. 

In order to generalize this result from spherical symmetry to 3-dimensions one additionally 
needs to show that the new uq, as a solution to the linear equation ffT2]) with the initial data 
( I13p . is near to the analogue of uq given by flHj) , containing now the Friedlander's radiation 
field (for the details see [g, Lemmas 2.3-2.5]). This step does not seem to essentially modify 
the rest of the proof except that spherical means will appear at various places and the 
constant A p will include integration over spheres, too. 

As already mentioned, the scaling technique developed here is not by itself restricted to 
the semilinear wave equations of the form ([I]). It can be easily generalized to any system 
of nonlinear wave equations. However, the strength of the estimate on the error term will 
crucially depend on the availability of corresponding pointwise decay estimates. 
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